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■ Abstract 

^ . The dispersion relation and mode amplitudes of obhque surface waves propagating on an acoustic 

O , 

tj^ ', double comb filter are obtained with a method based on the calculus of residues. We obtain a better 

O 

*S ' agreement (below 480 Hz) between theoretical predictions and measurements reported previously 

• when the filter was being supposed to be made of a single comb structure. 
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I. INTRODUCTION 



The behavior of a slow wave filter made of a pair of planar periodic waveguides subjected 
to low frequency acoustic waves incident upon the aperture separating the waveguides has 
been investigated theoretically and experimentally for its potential use in acoustic filtering 
devices Each waveguide has a comb structure consisting of a periodic array of blades 
perpendicular to a base plane (Figure 1). 

Using a mathematical model borrowed from the study of electrical filters, a filter having 
the same geometric structure of a single comb waveguide has been analyzed previously [1]. 
The dispersion relation, amplitude and phase as functions of frequency and wave number 
were derived and compared to experiment. In this work, we extend our previous theoretical 
results and consider the actual nature of the filter consisting of the two waveguides facing 
each other. We derive the dispersion relation and reflection (transmission) coefficients of 
surface waves propagating along any oblique wave number in the plane parallel to the comb 
structure base planes. 

Our calculations are based on a weak-coupling approximation and in the limit of 
large distance separating the two structures. This means the separation is much larger 
than the inter-blade distance. The blades are supposed to have a vanishingly small 
thickness and we neglect possible reflections from the planar base affecting the propa- 
gating modes, by direct analogy with the electromagnetic case [2]. This is equivalent to 
assuming a slot depth large with respect to the inverse lowest attenuation of the structure [2] . 

Our work is organized as follows: In section II, we discuss the geometry, propagating 
modes dispersion and amplitude relation for the surface waves. Section III covers the com- 
parison with the experimental results and the conclusion is in Section IV. 

II. DISPERSION RELATION, MODES AND AMPLITUDES 

Periodic arrays of slotted waveguides stacked to form a rectangular [3] or prismatic 
[4] structure are good candidates for reducing environmental noise (0.1 to 2 kHz). Their 



2 



properties have been analyzed theoretically and experimentally [1, 3, 4] such as their 
reflection scattering of sound waves harmful to the general population living near highways 
or other sources of damaging sources of low frequency noise. It is important to understand 
how these structures absorb, reflect, transmit or phase delay the incoming sound waves 
reaching them with arbitrary time dependent angles. For the rectangular structure, we 
have already undertaken such study from the experimental point of view as well as from the 
theoretical one. In this work, we set out to investigate a new type of structure introduced 
in detail in Ref. 1 theoretically and experimentally (Fig. 1). 

We have studied the dispersion relation of acoustic waves impinging on the structure at an 
arbitrary fixed angle in the base plane, and measured the sound refiection and transmission 
with respect to the incident angle. Our prior theoretical investigation took account of a 
single comb structure only. Here we extend it and deal with a symmetrical weakly- coupled 
double comb structure [5] in the limit ^ ^ 1 where b is half the distance between the tip 
of the blades belonging to each of the waveguides and d is the inter-blade distance in any 
waveguide (Fig. 1). 
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FIG. 1: Geometry of the double comb structure waveguide. 



Following our notation [1] , we write for the acoustic fields in region A (free space) keeping 
the symmetric modes only: 

CO 

^A{x,y,z)= J2 ^ne-^'^"^-^'^'cosh(Q;„y) (1) 

n=— oo 

where and r are the propagation constants along x and z and a,„ is the attenuation 
constant along y. The propagation constant (3o defining the fundamental mode is determined 
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from the propagation geometry (Fig. 1 of [6]). It is equal to where 9 is the angle, the 
surface wave vector makes with the x-axis [Fig. 1]. The surface wave has a smaller velocity 
than in true free space by the ratio V^5^~+T^/fc. In region B, the acoustic field in the n-th 
slot defined by the inequahties: ud — ^ < x < vd + ^ is given by: 

$^(x,y,^)= i?>-^'^^cos(— ^)cos[7„.(t/ + 6 + M] (2) 

n=-oo ^ 

The coefficients are determined with the help of Floquet's [7] theorem 5^ = 
BmG~^'^^°'^ and the abscissae are equal to x-(z/-l/2)d. In order to find the dispersion 
equation of the surface waves and the coefficients Am, B^, we will proceed as we did in our 
previous work following the approach pioneered by Whitehead [7] and Hurd [2]. It consists 
of writing the equations of continuity for the fields ^^and $b and their derivatives along the 
vertical y axis on the boundaries y — ±6. These equations are considered as originating from 
Cauchy's theorem of residues for a meromorphic function f{w) taken along some contour 
and the contribution of each pole is identified with the contribution of some corresponding 
mode. The contour and f{w) should be such that the presumed theorem of residues is 
satisfied. Moreover, the asymptotic behavior of f{w) is tailored by the underlying physical 
problem and is basically dictated by the scattering of the waves by the edges of the blades 
[7]. We obtain the following meromorphic function f{w) of the complex variable w: 

[eJ/^od_i] w-ao riiUTo) 112^ tt 
where Hil'"^) 3'ndn2('"^) the following infinite products: 

111 [w) = {[{w- j-fp){ )e (4) 

p=i P'^ 



and: 



U2iw) = Hiw - cyp)iw + a-p){- — feT^ (5) 
p=i 2p7r 



The propagation constants 7^ along y, are given by: 



7i = k'-T'-{^r withm=0,l... (6) 

In order to derive the dispersion relation, we form the ratio: 
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/(j7o) 

Taking the logarithm and using trigonometric identities [Ref. 2], we obtain: 



7o/i-7o = --sm — +2^ tan +- — + tan — 

71-2 /5o ^1 a-p 27rp |7p| pvr 

This equation is the same as that obtained by Hougardy and Hansen who treated 
a single comb structure from the electromagnetic point of view. Here, we are dealing 
with the weak coupling symmetric case limit and with the additional simplifying as- 
sumptions: ^ ^ l,aob ^ 1 and a_p ~ /5_p. We find that the dispersion relation is 
essentially the same as in the case of a single comb structure. The double comb structure 
simply behaves as a single one from the dispersion relation point of view. This justifies 
our assumptions in Ref. 1 where we found very good agreement between theory and 
experiment up to frequencies on the order of 400 Hz. Nevertheless this is not true for 
reflection (transmission) coefficients of the single/double comb structures as discussed below. 



In order to calculate the mode amplitudes and obtain from them the reflection (trans- 
mission) coefficients of the structure, we use the residue of f{w) at w = «„: 

Res[fiw)]^^^^ = A„/?„e^-^"'^/2 cosh(a„6) (9) 
to obtain (n=0, 1, 2...): 





1 ci7oe""^ 


a„ + ao\ 


a„ + ai 


1 (-^n 


- a^i\ 


r[2 + '^""] exp( ""'^ ^°(')) 


l^ol 


IGvr cosh(a„6) 






Ittn + j7l 


r[2 + ^(«„ + Po)] r[2 + ^(«„ - Po)] 



(10) 

where F stands for the Euler Gamma function. For negative values of n, it suffices to 
change a„ into — a„ in the above expression. Let us note that when the separation 2b 
between the two parts of the structure, becomes very large we recover exactly the expression 
found by Hougardy and Hansen |H] corresponding to a single comb structure. 
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In order to calculate the Br, coefficients, we use: 



^^"^ 2[(-)"exp(-j/3orf)-l] ^''^ 



and the definition (3) of fiw) to obtain: 



\Bn 




27oe 






1 ni(-J7.)l sin(f^) 


\Bo\ 






\hn + ao 




n2(-J7n)| (^) 



nnd ln(2) 
exp( ) (12) 

TT 



where e = 1 for n even, and e = hiZTrr- for n odd. 

|tan(^)l 



Let us note that the 5„ coefficients are the same as those obtained by Hougardy and 
Hansen P] refiecting the fact, the weak- coupling approximation affects in a different way 
the An and the Bn coefficients. This has important implications on our measurements of 
the amplitude profile. 

III. COMPARISON WITH EXPERIMENT 

In our previous work, we derived the dispersion relation, transmission and reflection 
coefficients and found excellent agreement between the single comb structure theory and 
experiment up to 400 Hz This work shows that a weak coupling between two comb 
structures does not affect the surface wave dispersion relations and the Bn amplitude 
coefficients but it does affect the An amplitude coefficients. 

We are going to evaluate how our theory modifies the amplitude ratio associated 
with the fundamental mode (n=0) in relation (10) compared with same given in Lahlou et 
al. [0. The double comb over single comb structure ratio of the two expressions is given by: 

For a given frequency and a given incident angle 6 we solve the dispersion relation given 
by equation (8), obtain the propagation factor ao and use it in (13). The corrections F(^) 
in dB are plotted versus 6 in the interval [1, 80] degrees for various frequencies [400-600 Hz] 
in Fig. 2. 
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FIG. 2: Corrections F{9) to the fundamental mode amplitude ratio with the following values 
(from the experimental setup) b=0.0125 m, d=0.05 m, h=0.112 m. The corrections calculated from 
10 logio{ ^p^1^ ^ evaluated as a function of the incident angle 9 at a fixed frequency varying 
from 400 to 600 Hz by steps of 50 Hz. 



The correction comprised between and -3 dB is small for higher frequencies and small 
incident angles. It decreases rapidly for angles larger than 10 to 20 degrees and by a larger 
amount for higher frequencies. A comparison to the experimental data reveals that the 
correction is pronounced mostly at higher frequency (336 Hz) and for the largest angle of 
incidence (47 degrees). For the highest experimental frequencies (480 and 496 Hz), the 
correction introduces more disagreement between the experimental and theoretical single 
comb structure theory. This behavior may be explained by the fact that there are several 
sources of errors associated with the measurements at these higher frequencies. 



IV. CONCLUSION 



We have developed a weak coupling theory based on the calculus of residues in order to 
model the oblique propagation of acoustic waves propagating through a slow wave filter made 
of a pair of comb structured waveguides separated by a distance that is large with respect to 
the inter-blade distance. The correction arising from the symmetrical coupling between the 
two waveguides has been evaluated and shown to improve slightly the agreement between 
the theoretical and the experimental values of Lahlou et al. being at the most 3 dB for 
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the largest frequency and angle evaluated. Those results show that the approximation taken 
in our previous investigation is quite acceptable and that the new theory does not bring 
substantial additional accuracy to our previous single comb structure model. Our studies of 

the strong coupling case (b < d) being mathematically much more complicated, and intended 
for improving the agreement between the theoretical results and the experimental ones at 
the higher frequencies are in progress and will be reported in the near future. 
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